ABSTRACT
Transmembrane protease serine 2 (TMPRSS2) has been identified as a critical key for the entry of coronaviruses into human cells by cleaving and activating the spike protein of SARS-CoV-2. To block the TMPRSS2 function, 18 approved drugs, containing the guanidine group were tested against TMPRSS2's ectodomain (7MEQ). Among these drugs, Famotidine, Argatroban, Guanadrel and Guanethidine strongly binds with TMPRSS2 S1 pocket with estimated Fullfitness energies of -1847.12, -1630.87, -1605.81 and -1600.52 kcal/mol, respectively. A significant number of non-covalent interactions such as hydrogen bonding, hydrophobic and electrostatic interactions were detected in protein-ligand complexes. In addition, the ADMET analysis revealed a perfect concurrence with the aptitude of these drugs to be developed as an anti-SARS-CoV-2 therapeutics. Further, MD simulation and binding free energy calculations were performed to evaluate the dynamic behavior and stability of protein-ligand complexes. The results obtained herein highlight the enhanced stability and good binding affinities of the Argatroban and Famotidine towards the target protein, hence might act as new scaffolds for TMPRSS2 inhibition. Communicated by Ramaswamy H. Sarma.
ABSTRACT
The SARS-CoV-2 virus is causing COVID-19 resulting in an ongoing pandemic with serious health, social, and economic implications. Much research is focused in repurposing or identifying new small molecules which may interact with viral or host-cell molecular targets. An important SARS-CoV-2 target is the main protease (Mpro), and the peptidomimetic α-ketoamides represent prototypical experimental inhibitors. The protease is characterised by the dimerization of two monomers each which contains the catalytic dyad defined by Cys145 and His41 residues (active site). Dimerization yields the functional homodimer. Here, our aim was to investigate small molecules, including lopinavir and ritonavir, α-ketoamide 13b, and ebselen, for their ability to interact with the Mpro. The sirtuin 1 agonist SRT1720 was also used in our analyses. Blind docking to each monomer individually indicated preferential binding of the ligands in the active site. Site-mapping of the dimeric protease indicated a highly reactive pocket in the dimerization region at the domain III apex. Blind docking consistently indicated a strong preference of ligand binding in domain III, away from the active site. Molecular dynamics simulations indicated that ligands docked both to the active site and in the dimerization region at the apex, formed relatively stable interactions. Overall, our findings do not obviate the superior potency with respect to inhibition of protease activity of covalently-linked inhibitors such as α-ketoamide 13b in the Mpro active site. Nevertheless, along with those from others, our findings highlight the importance of further characterisation of the Mpro active site and any potential allosteric sites.
Subject(s)
Antiviral Agents/pharmacology , Coronavirus 3C Proteases/antagonists & inhibitors , Coronavirus 3C Proteases/chemistry , Coronavirus Protease Inhibitors/pharmacology , Protein Multimerization/drug effects , SARS-CoV-2/drug effects , SARS-CoV-2/enzymology , Small Molecule Libraries/pharmacology , Amides/chemical synthesis , Amides/chemistry , Amides/pharmacology , Antiviral Agents/chemical synthesis , Antiviral Agents/chemistry , Azoles/chemical synthesis , Azoles/chemistry , Azoles/pharmacology , Coronavirus 3C Proteases/metabolism , Coronavirus Protease Inhibitors/chemical synthesis , Coronavirus Protease Inhibitors/chemistry , Humans , Isoindoles , Ligands , Lopinavir/chemical synthesis , Lopinavir/chemistry , Lopinavir/pharmacology , Microbial Sensitivity Tests , Models, Molecular , Molecular Structure , Organoselenium Compounds/chemical synthesis , Organoselenium Compounds/chemistry , Organoselenium Compounds/pharmacology , Ritonavir/chemical synthesis , Ritonavir/chemistry , Ritonavir/pharmacology , SARS-CoV-2/metabolism , Small Molecule Libraries/chemical synthesis , Small Molecule Libraries/chemistryABSTRACT
(1) Background: Virtual screening studies on the therapeutically relevant proteins of the severe acute respiratory syndrome Coronavirus 2 (SARS-CoV-2) require a detailed characterization of their druggable binding sites, and, more generally, a convenient pocket mapping represents a key step for structure-based in silico studies; (2) Methods: Along with a careful literature search on SARS-CoV-2 protein targets, the study presents a novel strategy for pocket mapping based on the combination of pocket (as performed by the well-known FPocket tool) and docking searches (as performed by PLANTS or AutoDock/Vina engines); such an approach is implemented by the Pockets 2.0 plug-in for the VEGA ZZ suite of programs; (3) Results: The literature analysis allowed the identification of 16 promising binding cavities within the SARS-CoV-2 proteins and the here proposed approach was able to recognize them showing performances clearly better than those reached by the sole pocket detection; and (4) Conclusions: Even though the presented strategy should require more extended validations, this proved successful in precisely characterizing a set of SARS-CoV-2 druggable binding pockets including both orthosteric and allosteric sites, which are clearly amenable for virtual screening campaigns and drug repurposing studies. All results generated by the study and the Pockets 2.0 plug-in are available for download.